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^ ; Abstract 

The equation of state of Hamiltonian lattice QCD at finite density is examined in the strong 
coupling limit by constructing a solution to the equation of motion corresponding to an effective 

o , 

I Hamiltonian describing the ground state of the many body system. This solution exactly diago- 

nalizes the Hamiltonian to second order in field operators for all densities and is used to evaluate 
^ . the vacuum energy density from which we obtain the equation of state. We find that up to and 



beyond the chiral symmetry restoration density the pressure of the quark Fermi sea can be negative 
^ ' indicating its mechanical instability. Our result is in qualitative agreement with continuum models 

and should be verifiable by future lattice simulations of strongly coupled QCD at finite density. 
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I. INTRODUCTION 



Lattice gauge theory is currently the only known method of solving Quantum Chromo- 
dynamics (QCD) from first principles. It has developed sufficiently to be able to calculate a 
broad range of low and intermediate energy hadronic observables from ground state hadron 
masses to pion-nucleon scattering lengths. In addition, lattice studies of QCD at finite tem- 
perature (T), especially its spectacular success in demonstrating the deconfinement phase 
transition, have been invaluable and continue to play an important role in the search for the 
quark-gluon plasma. 

However, as is well known, one of the outstanding problems in lattice gauge theory is 
the consistent implimentation of chemical potential in numerical simulations [jll. Progress 
in lattice QCD calculations at finite chemical potential (yu) with dynamical quarks has 
been hindered by the presence of the complex fermion determinant which renders standard 
Monte Carlo techniques useless. In fact, currently there is only one numerical method of 
simulating finite density QCD with three colors at zero temperature. This method is based 
on the Monomer-Dimer-Polymer algorithm developed by Karsch and Miitter However 
its applicability is limited to the strong coupling limit, and furthermore a recent study 
indicates that this algorithm might not be reliable for studying the chiral phase transition 
at finite density. Therefore even a qualitative description of finite density lattice QCD is 
welcome. 

One method of studying finite density QCD on the lattice is to invoke the strong cou- 
pling approximation where analytical methods are applicable. Although far from the re- 
alistic continuum limit, the strong coupling approximation has played an important role 
in the development of QCD lattice gauge theory from its very inception. In the renowned 
paper by Wilson [^] this approximation was invoked to demonstrate quark confinement on 
the Euclidean space-time lattice. Soon thereafter Kogut and Susskind § formulated the 
Hamiltonian lattice gauge theory and concluded that in the strong coupling limit the quark 
dynamics is best described by a collection of non-Abelian electric flux tubes with quarks 
attached at their ends. This was followed by the work of Baluni and Willemsen who used 
a variant of the Kogut-Susskind formalism to demonstrate quantitatively that dynamical 
chiral symmetry breaking indeed takes places in lattice QCD at strong coupling [Q. Finally, 
calculations by Kogut, Pearson and Shigemitsu [0] and by Creutz [|] suggesting the absence 
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of a phase transition between the strong and weak couphng regimes of QCD motivated 
numerous studies using the strong couphng approximation. 

Strong couphng QCD at finite T and/or /i has previously been studied analyticaUy both 
in the Euchdean ||^, |TD|, |rT|] and in the Hamihonian |]T^, |T^, |T^, |T3| formulations. One of the 



main objectives of these studies was to investigate the nature of chiral phase transition at 
finite temperature and density. In each study this was accomplished by constructing some 
effective action or Hamiltonian for strongly coupled lattice QCD using Kogut-Susskind 
fermions. These effective descriptions involve the introduction of composite meson and 
baryon fields which are treated in the mean field approximation.^ The consensus is that at 
zero or low temperatures strong coupling QCD undergoes a first order chiral phase transi- 
tion from the broken symmetry phase below a critical chemical potential fic to a chirally 



symmetric phase above ^c- The only exception is the work by Le Yaouanc et al. |14| which 
does not involve effective composite fields but is equivalent to the mean field approximation. 
In this case the chiral phase transition was found to be of second order. 

In this paper we present the equation of state of Hamiltonian lattice QCD at finite density 
in the strong coupling limit using both Kogut-Susskind and Wilson fermions. As in previous 
studies we begin with an effective theory by using a Hamiltonian describing the ground state 
of strongly coupled QCD. However, our approach differs from earlier works in that we do 
not introduce composite fields but explictly construct a solution to the field equations of 
motion corresponding to the effective Hamiltonian. This solution exactly diagonalizes the 
Hamiltonian to second order in field operators for all densities and is used to calculate the 
vacuum energy density from which we obtain the equation of state. We find that up to 
and beyond the chiral symmetry restoration density the quark Fermi sea can have negative 
pressure indicating its mechanical instability. Our result is in qualitative agreement with 



those obtained using continuum effective QCD models [|T6|, and should be verifiable by 
future lattice simulationns of strongly coupled QCD at finite density. 

Our approach admits to a systematic extension to finite temperature and to the descrip- 
tion of bound states. In fact, we first introduce temperature and chemical potential simul- 
taneously into our formalism and then take the limit of vanishing temperature to examine 
the consequences. Description of bound states is accomplished by interpreting our solution 



^ The Monomer-Dimer-Polymcr algorithm M also uses composite meson and baryons fields. 
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within the context of the N-quantum approach (NQA) to quantum field theory [|T8], [19 



which we shall discuss in the concluding section. In the same section we propose how the 
NQA may be combined with the present approach to study the nature of the deconfinement 
phase transition. 

In Section^ we introduce our effective Hamiltonian for the ground state of strong coupling 
QCD using Wilson fermions, and discuss the condition under which it can be extended to 
finite T and fi. Our ansatz for the lattice quark field at finite T and /i is presented in 



Section |T|. The equation of motion at finite /i is then derived in Section |I^ and used to 
diagonalize the effective Hamiltonian to second order in field operators and to evaluate the 
vacuum energy density. In the same section we determine the unknown quantities in our 
ansatz by deriving coupled equations for the dynamical quark mass and the total chemical 
potential and solving them self-consistently. Having constructed a solution for the quark 
field we present in Section the equation of state of Hamiltonian lattice QCD at finite 
density in the strong coupling approximation. We summarize our results in Section ^ and 
discuss how our approach may be extended to incorporate temperature and to describe the 
deconfinement phase transition. A review of the properties of free lattice Wilson fermions 
using the Hamiltonian formulation is given in the Appendix. 

II. THE EFFECTIVE HAMILTONIAN 



We begin by introducing Smit's effective Hamiltonian ||2^ describing the ground state of 
strongly coupled QCD. This state is the one in which no links are excited by the electric 
flux. It is also inflnitely degenerate since various color singlet states may be created at each 
lattice site without increasing the ground state energy. This degeneracy is lifted by the 
propagation of quarks on the lattice. The simplest type of such a propagation involves a 
quark exciting a flux link and an anti-quark deexciting the same link and corresponds to the 
propagation of a meson. Smit obtained an effective Hamiltonian describing this propagation 
using second order perturbation theory involving only the quark fleld \1/ with a nearest 
neighbour interaction. The Hamiltonian is effective because it only acts on the space of 
states with no excited links. Nevertheless, it serves our purpose since the main quantity of 
interest in this work is the vacuum energy density which is obtained by diagonalizing Smit's 
Hamiltonian. 
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In the Hamiltonian formulation of lattice field theory only the spatial coordinates are 
discreticized while the temporal coordinate remain continuous. We adopt the notation of 
Smit 1^ where the discrete sums over the spatial and momentum coordinates are given by 



E^^'E (1) 

X fa 

and 

e4? 

p n 

with X = am = a(mi, 7712,7713) and p = j^n = -^(71-1,772,773). Here a is the lattice spacing 
and L defines the number of unit lattice cells with mi,ni = 0, ±1, ±2, ±L. With this 
notation the volume V is given by = (2La)^. Henceforth we shall work in lattice units 
where a = 1 so that — tt < pi < +tt, and use Wilson fermions. 



The effective Hamiltonian derived by Smit [20| using the temporal gauge is 



= Mo5^(vI/L)^(f)(^o)^^(vl/,,)^(f) 



X 

—EE 

X I 



+ (vl>y (f +aO (EO;, {'^ap)six){^L)(x) (E,),, {^ba),ix + d,) 



(3) 



where E/ = —7(707; — ir'-fo) and di is a unit vector along the positive /-axis. We denote 
color, flavor and Dirac indices by (ab), (ajS) and {pwyS), respectively. Summation convention 
for repeated indices is implied. The three parameters in this Hamiltonian are the Wilson 
parameter r which takes on values between and 1, the current quark mass Mq and the 
effective coupling constant K = 2Nc/{N^ — l)l/g^ where g is the QCD coupling constant. 
Nc is the number of colors. When r = the quark fields become Kogut-Susskind fermions. 

Smit's Hamiltonian is valid to order Oll/g"^) in the strong coupling expansion. The 
0(1/(7^) corrections involve products of quark bilinears which describe meson propagation 
mentioned above and are known as "meson terms". For A^c = 3, contributions from the 
subsequent order in the l/g"^ expansion would consist of products of terms which are trilinear 
in the quark fields called "baryon terms". These meson and baryon terms appear in the 



strong coupling expansions of both Euclidean and Hamiltonian lattice QCD and are the 
motivations for introducing effective composite meson and baryon fields. In this work we do 
not take the baryon terms into account but our formalism presented here is also applicable 
if such terms were present in the effective Hamiltonian. 

In the absence of the current quark mass and the Wilson parameter (Mq = r = 0), 
Hes posseses a U{4:N{) symmetry with A^f being the number of flavors. This symmetry is 
spontaneously broken to U {2Nf)®U{2Nf) accompanied by the appearance of 8Nf Goldstone 
bosons A finite current quark mass also breaks the original U{4Nf) symmetry, albeit 
explicitly, to U{2Nf) ® U(2Nf). Introduction of the Wilson term explicitly breaks the latter 
symmetry further down to U{Nf) thereby solving the fermion doubling problem. 

We shall work exclusively in momentum space. Our convention for the Fourier transform 
from configuration to momentum space is \E'(x) = Yl^'^ip)^^^ ^ which implies that the 
volume V is given hj V = = ^p-^p-. Then the charge conjugation symmetric form of 
Smit's Hamiltonian in momentum space is given by 



This effective Hamiltonian is the starting point of the present investigation. Our method 
for obtaining the equation of state consists of extending H^s to finite n and constructing a 
quark field operator \E' which diagonalizes the Hamiltonian to second order in field operators 
for all densities. Once this solution has been found it can be used to evaluate the vacuum 
energy density from which we obtain the pressure of the many body system. 

However before extending H^s to finite T and/or /i, it is necessary to impose a condition 
on these external parameters so that all links would remain in their ground states. In the 
strong coupling limit the amount of energy required to excite one color electric flux link is 



-5^Mo(7o),. (*L),(P),(M.(-P) 



p 




(4) 




K 



1 



(5) 
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Therefore an extension of Hcs to finite T and/or fi will be valid as long as T, fi < 1/K [p^^ 
since the Hamiltonian only acts on the space of states with no excited links. We shall see 
that this condition is satisfied in the present work. 

The effective Hamiltonian is extended to finite T and fi in two steps. The first one is to 
make the following trivial replacement of the current quark mass term in H^s 

Mo(7o)p,-^Mo(7o)p,-/io5p. (6) 

where /io is the quark chemical potential. Note that fio should not be identified with the 
total chemical potential //tot of the interacting many body system. As we shall see below, the 
O {1/ g"^) interaction terms in Hes will induce a correction to fiQ which in general is momentum 
dependent. We shall therefore refer to /io as the "bare" quark chemical potential and treat 
it as an input parameter. The second step is to introduce an ansatz for the quark field at 
finite T and fi. 



III. AN ANSATZ FOR FINITE T AND ^ 

We proceed by presenting our ansatz for the \E' field in H^s at any temperature and 



chemical potential. The special case of this ansatz for free space was given in ||2T|. It has the 
same structure as the free lattice Dirac field and obeys the free lattice Dirac equation with 
a mass which is interpreted as the dynamical quark mass. This mass is the only unknown 
quantity in the free space ansatz and is determined by solving a gap equation. It was shown 
in that this ansatz exactly diagonalizes Hcs to second order in field operators. Properties 



of free lattice Dirac fields using Wilson fermions are given in the Appendix. 



Temporarily dropping color and flavor indices the free space ansatz given in [21 



IS 



with V denoting the Dirac index. The annihilation operators for particles h and anti-particles 
d annihilate an interacting vacuum state \Qq), and obey the free fermion anti-commutation 
relations. The properties of the lattice spinors ^ and r] are given in the Appendix. The free 



Note that in Wm E has been approximated hy E k 
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lattice Dirac equation fixes the excitation energy uj{p) to be 



1/2 



uj{p) = i ^sin^(p- ni) + M'^{p) 



(8) 



where M{p) is the dynamical quark mass. 

In order to extend Eq. (0) to finite T and fi we observe that the annihilation operators b 
and d in v]/^''^ no longer annihilate the interacting vacuum state at finite T and fi denoted as 
\Q(T,fi)). To construct operators that annihilate | Q(T,^)) we apply a generalized thermal 
Bogoliubov transformation to the b and d operators following the formalism of thermal field 
dynamics 



bip) = a,B{p) - f3,B\-p) 
d{p) = ^,D{p) - 5,D\-p) 



(9a) 
(9b) 



The thermal field operators B and W annihilate a quasi-particle and create a quasi-hole at 
finite T and /i, respectively, while D and D'^ are the annihilation operator for a quasi-anti- 
particle and creation opertor for a quasi-anti-hole, respectively. 

These thermal annihilation operators annihilate the interacting thermal vacuum state for 
each T and fi. 



B{p)\giT,fi)) = Bip)\giT,fx)) = Dip)\g{T,fx)) = Dip)\giT,fi)) =0 



(10) 



We note that the thermal doubling of the Hilbert space accompanying the thermal Bogoli- 
ubov transformation is implicit in Eq. ( |T0|) where a vacuum state which is annihilated by 
operators B, B, D and D is defined. In addition, since we shall be working only in the 
space of quantum field operators it is not necessary to specify the structure of the thermal 
vacuum \g(T, ft)). 

The thermal operators also satisfy the Fermion anti-commutation relations 



B\p),B{q) 



D\p),D{q\ 



B\p),B{q 



D\p),b{q 



:ii) 



with vanishing anti-commutators for the remaining combinations. The coefficients of the 
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transformation are 



1 — nz 



(12a) 
(12b) 
(12c) 
(12d) 



where 



n: 



± 



(13) 



"P ^{uj(p)±t,)/T _^ I 

are the Fermi distribution functions for particles (n~) and anti-particles (^^ ). We stress that 
the chemical potential appearing in the Fermi distribution functions is the total chemical 
potential of the interacting many body system. The coefficients are chosen so that are 
given by 



n- = {g{T,f,)\b\p)b{p)\g{T,^^)) 



n; = {g{T,^^)\d\p)d{p)\Q{T,^i)) 



(14a) 
(14b) 



Hence in this approach temperature and chemical potential are introduced simultaneously 
through the coefficients of the thermal Bogoliubov transformation and are treated on an 
equal footing. 

After applying the Bogoliubov transformation to Eq. (|^) our ansatz at finite T and n 
becomes 



^.(t,p) 



a. 



+ 



^pD\-p) - SpD{p) 



(15) 



and satisfies the equation of motion corresponding to the free lattice Dirac Hamiltonian at 
finite chemical potential given by 



- ^ sin(p ■ ni){'yo'ji)nu + M{p){'yo)nu - yUtotV 



^l{t,p),^,{t,p: 



(16) 



The spinors ^ and t] in Eq. (|15D obey the same properties as in free space and the excitation 
energy uj{p) has the same form as in Eq. (H). The unknown quantities in our ansatz Eq. (plSf ) 
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/^tot 



are the dynamical quark mass M{p) and the total chemical potential ^tot which will be 
determined in the following section. 

In this work we shall take the T — limit which amounts to setting 7p = 1 and 6p = 
in the Bogoliubov transformation Eq. (^) thereby suppressing the excitation of anti-holes. 
In this limit becomes the Heaviside function = 9{fitot — ^ij>)) defining the Fermi 
momentum pp through the relation 

\ 1/2 

Y,s\ri\pF-n{) + M\pF)\ (17) 

Note that we define chemical potential such that /i^ot > M{pp) which differs from the 
conventional definition of chemical potential used in lattice calculations where > 0. 

One of the simplest quantities to calculate using the ansatz of Eq. ( [I5| ) in the T — > 
limit is the quark number density n given by 

c ^ p 

p p 

where the symbol : : denotes normal ordering with respect to the vacuum at zero temper- 
ature I Q(T = 0,yu)). Therefore, above a sufficiently large value of /itot the quark number 
density becomes a constant which with the present normalization will equal unity. This 
saturation effect is purely a lattice artifact originating from the sin^(p ■ hi) term in uj{p). 

Another quantity that may be readily calculated using the T — > ansatz is the chiral 
condensate. It is found to be proprotional to the dynamical quark mass 

Below we shall derive a gap equation for M{p) and show that for a given physically rea- 
sonable set of parameters there exists a critical chemical potential above which M{p ) = 0. 
Thus the chiral condensate may be identified as being the order parameter for the chiral 
phase transition at finite density. 
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IV. APPLICATIONS OF THE EQUATION OF MOTION 



A. The equation of motion 

We now calculate the equation of motion corresponding to i^cfr with our ansatz for finite 
fi using two light flavors. The result is used to show that our ansatz exactly diagonalizes the 
Hamiltonian to second order in field operators for all densities and to calculate the vacuum 
energy density. In addition, by analyzing the Dirac structure of the equation of motion we 
derive coupled equations for the dynamical quark mass and the total chemical potential. 
They are solved to lowest order in the l/N^ expansion thereby completing our construction 
of a solution to the lattice field theory defined by H^s- 

The equation of motion for Hes is obtained by exploiting the fact that our ansatz also 
satisfies the equation of motion corresponding to the free lattice Dirac Hamiltonian given 
in Eq. (|T^) . We therefore have the relation 

(20) 

which plays a crucial role in our construction of a solution for the quark field Evaluating 
both sides of Eq. (|20|) and equating terms which are linear in the field operators we obtain 
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the equation of motion for \I/ 
I 

{P- q)-ni ({^i)J^i)ls + (sol. (SO, 5] 
+ COS {p + q )-h, (^(S0;(S0I5 + (S0p.(s0,5 



7(5 



+ COS {p + q).h,\ + 



{^ao),{t,q) (21) 



with A+(p) = S,{p) ® i'^ip) being the positive energy projection operator defined in 
Eq. ( lATSal) . 



B. Diagonalization of i/eff and the vacuum energy density 

We shaU now show that our T ^ ansatz exactly diagonahzes the effective Hamihonian 
to second order in field operators. The diagonalization procedure involves only algebraic 
substitutions and does not require any approximations. The quantity of interest here is the 
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off-diagonal Hamiltonian wliicli, to second order in field operators, is found to be 



+ cos ip + q)-ni(^{^i)l{^i)ls+i^i)J^d,^ 



1 K 

XT 



® eKg ) 



+ COS (p + g ) ■ n,i^{j:,)l{j:,)l^ + 



'Vs[-q. 



(22) 



We see from Eq. (|2^) that the elementary excitations of the effective Hamiltonian are color 
singlet (quasi) quark-anti-quark excitations coupled to zero total three momentum. They 
correspond to the meson propagation on the lattice responsible for lifting the degeneracy of 
the ground state of strongly coupled QCD. 

With the use of the equation of motion for the \l/ field Eq. ( pT] ) , the equation of motion for 
the rj spinor Eq. ( [Al(j| ) and the orthonormality condition for the ^ and rj spinors Eq. 
we can show that 



Hos\g{0,fi)) = ^<ag^J(g)[- J^sin(g-fi,)(7o70p5-M(g)(7o)p5 + /itot5p5 

q I I 

= El l^i'l) + /^totl r/p(-g) \Bl^{q)Dl^{-q)\g{0,fi)) 



vs[-q. 







(23) 



Note that this result is valid for any dynamical quark mass and total chemical potential. 
Therefore in the T ^ limit our ansatz shown in Eq. (|l^) exactly diagonalizes the effective 
Hamiltonian to second order in field operators for all densities. 
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Having diagonalized the second order Hamiltonian we can proceed to evaluate the vacuum 
energy density. Using Eq. (|21| ) once more we find 



a, 



;MoTr(A+(p)7o)+2/3pVo 



p,q I 



+ cos ip + q).J{^,)l{i:,)\^+{j:,)Jj:,)^^ 



K 



p,q I 



+ cos ip + q).n,{ + 



p,i 



2 uj{p) 



(1 - r^) cos(p + q)-hi- fitot 



+ (l + /?>o (24) 



For free space the diff'erence of the vacuum energy densities in the Wigner-Weyl (M(g ) = 0) 
and Nambu-Goldstone (M(g ) 7^ 0) phases of the theory is positive 



(25) 



Numerically we find that Eq. (^) also holds for finite chemical potential. Therefore the 
true ground state of our interacting many body system is in the phase with broken chiral 
symmetry. 



C. Dynamical quark mass and ^tot 

We now derive the equations for the dynamical quark mass and the total chemical poten- 
tial and solve them to determine our solution Eq. (|l^) for each density at zero temperature. 
To accomplish this we explicitly evaluate the right hand side of Eq. (^1]) to reveal its Dirac 
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structure. The result may be cast in the following compact form 

J^sin(g ■ hi){-foji)^s + M(g)(7o)^5 - fJ'tot^uS (^aaj^l^,?) = 
I 

A{q){^oli).5 + B{q){^o)uS + C{q)6,s\ {^aa)s{t, q) (26) 

The equations for M{p) and /itot are obtained by equating the coefficents of the 70 operator 
and the Kronecker delta function, respectively. 

The gap equation determining M{p) is given by the coefficient B{q) 

M{q) = B{q) 



Mo+^K(i-ox:(i-/5,^) 

Ir'^ cos(p ■ hi) cos(g ■ hi) — -(1 + r^) cos(p + q) ■ hi^ (27) 

The structure of this gap equation is very similar to the one in free space {(3^ = 0) found in 
|2l| . The dynamical quark mass is a constant to lowest order in Nc but becomes momentum 



dependent once 1/Nc correction is taken into account. 

Similarly, the total chemical potential is given by the coefficient C{q) 

A^tot = -C{q) 

= -"0 + E E K + - 2 (1 - r') cos {p + q)] (28) 

Thus /itot is a sum of the bare chemical potential fiQ and an interaction induced chemical 
potential which is proportional to the effective coupling constant K. Furthermore, the latter 
contribution to yUtot is momentum dependent and this dependence is a l/N^ correction just 
as in the case of the gap equation. It should be noted that the above shifting of the bare 
chemical potential by the interaction is not a new effect. For example, in the well-known 
and well-studied Nambu-Jona-Lasinio model pSf at finite T and /i the interaction induces 



a contribution to the total chemical potential which is proportional to the number density 



24, 25 



The two equations Eqs. (^7^ and (^8]) are coupled and therefore solutions for M and /itot 
must be found self-consistently for each value of the input parameter fj,Q. We shall solve the 
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FIG. 1: Dynamical quark masses in free space M(Free) to 0{N^) as functions of the effective 
coupling constant K obtained using Wilson parameters r = 0.00, 0.50 and 0.75. Critical coupling 
constants are Kq = 0.732, 0.976 and 1.673 for r = 0.00, 0.50 and 0.75, respectively. 

coupled equations to 0{Nl!) which is equivalent to invoking the mean field approximation. 
At this order in both the dynamical mass and the total chemical potential are momentum 
independent. It is also the same order in the 1/Nc expansion used to obtain results in all 
previous studies of strongly coupled lattice QCD. All results are presented using Mq = 
and Nc = 3. 

We first discuss the solutions to the gap equation in free space. In Figure |l] we show 
dynamical quark masses in free space M(Free) as functions of the coupling constant K for 
Wilson parameters r = 0.00, 0.50 and 0.75. The dynamical quark masses were obtained in 
a straightfoward manner by solving the free space gap equation 



M(Free) = Mq + -K{1 



M(Free) 
uj{p) 



(29) 



The figure shows that for each value of r there exists a critical coupling constant Kq > 
above which the theory is in the broken symmetry phase. This is also true for the r = case 
corresponding to the use of Kogut-Susskind fermions. In this case the symmetry breaking 
takes place only for K > 0.732. 

The dependence of the dynamical mass, and consequently of the chiral condensate through 
Eq. (|l^), on the coupling constant is quahtatively different from the results obtained previ- 



ously using the same effective Hamiltonian pO, p6[. In both |20| and pq] qq pair condensation 
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FIG. 2: Dynamical quark masses M as functions of the bare chemical potential to lowest 
order in N^. The dynamical mass labelled "Uncoupled" was obtained by simply solving Eq. (|2l 
with /itot = Mo and exhibits a first order phase transition with a critical chemical potential of 
(mo)c = 0.825. The result labelled "Coupled" was obtained by solving the coupled equations 
Eqs. ( p7| ) and (|2^) self-consistently. A first order phase transition also takes place, but now the 
value of (/Uo)c is 0.785. The Wilson parameter and the coupling constant are set to r = 0.0 and 
K = 0.9, respectively. 

occurs for any value of K > 0. We find that the attraction between a quark and an anti- 
quark must be sufficiently large enough for a qq pair to condensate in the vacuum. Thus 
our approach provides a mechanism for chiral symmetry breaking which other approaches 
do not. In addition, our results are consistent with the works by Finger and Mandula ||2^ 
and by Amer, Le Yaouanc, Oliver, Pene and Raynal ||28| who have shown that in QCD in 
the Coulomb gauge qq condensation takes place only above a critical coupling constant. 

Examples of finite /i solutions to the coupled equations Eqs. ([27| ) and ( p8l) to lowest order 
in Nc are shown in Figures ^ and |^. In Figure ^ we show dynamical masses as functions of 
the hare chemical potential /io to highlight the importance of solving the coupled equations 



consistently. The figure shows the dynamical quark mass obtained by solving only Eqs. (^Tf) 
with /itot = /^o as well the mass obtained by solving the coupled equations consistently. 
Using r = and K = 0.9 a first order phase transition is observed in both cases, but the 
values of the critical /zq are 0.825 when //tot = A^o and 0.785 when the two equations are 
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FIG. 3: Dynamical quark mass M as a function of total chemical potential ^Utot for two values of 
the effective coupling constant K. These results were obtained by solving Eqs. ( |27| ) and ( p^ ) self- 
consistently to lowest order in Nc using r = 0.25. There is a second order chiral phase transition 
when the effective coupling constant K is 0.9 with a critical chemical potential of (^tot)c ~ 0.716. 
The order of the phase transition becomes first order with (/itot)c ~ 0.871 when K is increased to 
1.0. 

solved self-consistently. Therefore the critical chemical potential will be overestimated if 
interaction induced corrections to the bare chemical potential are ignored. 

In Figure ^ we present the dynamical mass as a function of the total chemical potential 
/itot for two values of K obtained with r = 0.25. From the figure we see that the phase 
transition can be either first or second order depending on the value of the coupling constant. 
When K = 0.9 we find a second order phase transition with a critical chemical potential of 
(/Utot)c ~ 0.716, while if the coupling constant is increased to = 1.0 the phase transition 
becomes first order with a larger critical chemical potential of (/itot)c ~ 0.871. This increase 
in the critical chemical potential with K has also been observed in |]lT|. Furthermore, we 
find that when K = 0.9 lattice saturation sets in above (/itot)c at around /itot ~ 0.898 while 
this effect takes place immediately above (/Utot)c for K = 1.0. These values of chemical 
potentials are smaller than the energy E = 1/K required to excite one color electric flux 
link as given in Eq. (^. Therefore with a reasonable set of parameters it is possible to 



extend Smit's effective Hamiltonian to finite density as was first pointed out in [14 . 
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FIG. 4: Quark number density n as a function of total chemical potential ^tot for two values of 
effective coupling constant K with r = 0.25. When K = 0.9 there is a jump in the number density 
at the phase transition point at (/itot)c ~ 0.716 from n ~ 0.037 to n ~ 0.070, while for K = 1.0 n 
becomes unity immediately above the critical chemical potential of (//tot)c ~ 0.871 due to lattice 
saturation. 

Having solved the self-consistency equations for the dynamical quark mass and the total 
chemical potential to lowest order in A^c we have constructed a mean field solution for the 
quark field appearing in the effective Hamiltonian Eq. (1). In Figure | we show the quark 
number density obtained with this solution as a function of fitot for K = 0.9 and 1.0. In 
both cases the number density is a monotonically increasing function of yUtot in the broken 
symmetry phase. When K = 0.9 there is a jump in the number density at the phase 
transition point at (/itot)c ~ 0.716 from n ^ 0.037 to n 0.070. Beyond this point the 
number density continues to increase monotonically until when the lattice saturation sets in 
at /itot ~ 0.898. This behaviour of the number density is qualitatively the same as the one 
obtained numerically using the Monomer-Dimer-Polymer algorithm as can be seen from 
a comparison with Figure 5 of 0. For K = 1.0 the lattice saturation takes place at the 
phase transition point at (/itot)c ~ 0.871 and beyond this point the number density remains 
a constant at n = 1. Noting that the number density at the phase transition point is 
n ^ 0.013, the number density for = 1.0 may be approximated by a Heaviside function 
of the form n = 6{^to\. — (/^tot)c)- This is exactly the result obtained in Eq. (2.47) of ||T5[ 
where a different effective Hamiltonian was used to study strongly coupled lattice QCD 
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FIG. 5: Pressure as a function of total chemical potential fitot obtained using Kogut-Susskind 
fermions (r = 0.0) with K = 0.8 and 0.9. 




r = 0.00 



at finite density. Furthermore the functional form of the number density found in is 



independent of the strength of the interaction. Therefore the results presented in for 
> /^c represent nothing but those obtained in the lattice saturation limit. 



V. EQUATION OF STATE 

We are now in a position to determine the equation of state by numerically evaluating 
the thermodynamic potential density using the mean field solution determined above in the 
vacuum energy density Eq. (|2^) . In Figure ^ we plot pressure as a function of /itot for K 
= 0.8 and 0.9. The value of the Wilson parameter is r = 0.0 so that the results have been 
obtained using Kogut-Susskind fermions. For both values of K we find that the pressure 
of the quark Fermi sea is negative and monotonically decreasing in the broken symmetry 
phase. For = 0.8 the pressure remains negative but increases in the symmetry restored 
phase, at least until the lattice saturation point, and has a cusp where the two phases meet. 
Unfortunately, for K = 0.9 we can not make a definite quantitative statement concerning 
the behaviour of the pressure in the symmetry restored phase due to lattice saturation, 
except to mention that there is a discontinuity when going from one phase to another. We 
find qualitatively similar results when Wilson fermions are used to calculate the pressure as 
shown in Figure ^ The parameter used in this figure are r = 0.25 and K = 0.9 and 1.0. We 
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FIG. 6: Pressure as a function of total chemical potential /xtot obtained using Wilson fermions (r 
= 0.25) with K = 0.9 and 1.0. 

may therefore conclude, at least at the mean field level, that up to and beyond the chiral 
symmetry restoration point the quark Fermi sea can have negative pressure and therefore 
can be mechanically unstable with an imaginary speed of sound. 

Our conclusion regarding the (strongly coupled) quark matter stability at finite density 
is consistent with similar studies using effective continuum models of QCD. In the Nambu- 
Jona-Lasinio model |T6tl and the instanton induced 't Hooft interaction model [ITTi, mean 



field calculations show that cold and dense quark matter may be unstable in the phase with 
spontaneously broken chiral symmetry, but can become stable in the symmetry restored 



phase at high enough density. In particular, the result for the pressure obtained in [jT^ is 
qualitatively the same as the one shown in Figures ^ and ^ as can be seen by comparing 
the figures with Figure 1 of [|1^ . The possibility of unstable quark mattter lead the authors 



of [|T^ and [|T^] to speculate the formation of nucleon droplets, reminiscent of the MIT bag 
model, in the broken symmetry phase. We shall not indulge on such a speculation here since 
we are working in an artificial strong coupling regime. Nevertheless, our results concerning 
the negative pressure is certainly verifiable in future lattice simulations of finite density QCD 
at strong coupling. 
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VI. CONCLUSION AND OUTLOOK 



In this work we studied the equation of state of two flavored Hamiltonian lattice QCD in 
the strong coupling limit at finite density using both Kogut-Susskind and Wilson fermions. 
Starting from an effective lattice Hamiltonian for the ground state of the strongly coupled 
QCD, we constructed a mean field solution which exactly diagonalizes the Hamiltonian to 
second order in field operators for all densities. This solution obeys the free lattice Dirac 
equation with a dynamical quark mass and total chemical potential which are determined 
by solving a coupled set of equations obtained from the equation of motion. From the gap 
equation determining the dynamical quark mass we find that at the mean field level the 
order of the chiral phase transition can be either first or second order depending on the 
values of input parameters. 

The equation of state was obtained by evaluating the thermodynamic potential density 
from the vacuum energy density using our solution. We find that the pressure of the strongly 
interacting many body system may be negative in the broken symmetry phase indicating 
the mechanical instability of our quark Fermi sea. There are indications of this instability 
beyond the phase transition point although no definite conclusions could be reached for 
very high densities due to lattice saturation. Nevertheless this behaviour of the pressure 
was found both for the case of Kogut-Susskind and Wilson fermions and seems, at least at 
the mean field level, to be robust. In addition, our result concerning negative pressure is 
in qualitative agreement with studies using continuum effective QCD models, and therefore 
should certainly be verified by future lattice simulations of strongly coupled QCD at finite 
density. 

To include temperature into our formalism we simply repeat our calculations using the 
ansatz given in Eq. ([15D at non-zero T. Preliminary calculations indicated that, in addition 
to particle-anti-particle excitations, the elementary excitations would now involve particle- 
hole, anti-particle-anti-hole and hole-anti-hole excitations. Because of these additional 
types of excitations our ansatz would no longer be able to exactly diagonalize the second 
order Hamiltonian. In fact, a simple exercise would show that at finite T even the free lattice 
Dirac Hamiltonian Eq. (p!6D is not diagonal due to particle-hole and anti-particle-anti-hole 
excitations. 

We now turn our attention to the possibility of studying the nature of the confinement- 
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deconfinement phase transition. Our solution presented in this work is non-confining and 
therefore it would be hopeless to use it to study this important phase transition. What is 
lacking in our formalism is the description of bound states. However, our solution presented 
here is by no means unique or complete and it can be systematically improved to include all 
the bound states allowed by the effective Hamiltonian. This is accomplished by interpreting 
our solution within the context of the N-quantum approach (NQA) to quantum field theory 

MM- 

NQA is a method to solve field equations of motion by expanding the interacting Heisen- 
berg fields in terms of asymptotic fields obeying the free field equations of motion. Here the 
on-shell masses can but need not equal the physical masses of the fields. This expansion is 
known as the Haag expansion p9| and our ansatz presented here is nothing but the first term 
in this expansion. Note that because we are working in the Hamiltonian (Kogut-Susskind) 
formulation of lattice field theory the time variable is continuous and therefore we can in- 
troduce and work with the concept of asymptotic fields. The second order terms in the 
Haag expansion would consist of a product of fermionic quark fields and bosonic elementary 
color singlet qq bound state fields. The coefficient of each of the second order terms are 
interpreted as creation amplitudes for the bound states and are known as Haag amplitudes. 

Supressing color and fiavor indices for simplicity, our extended ansatz for the Heisenberg 
quark field \1/ to second order in the Haag expansion in free space will have the following 
structure 

^Free^^^ = ^l{q) + Y. [ S kd% 5\q + k - h) f^^ik ,h) : vl/J(-fc)4)(6) : (30) 

where the symbol : : denotes normal ordering. In Eq. (pO]), B is the elementary bosonic 
field while the superscript indicates that the fields obey their corresponding free field 
equations of motion. The Haag amplitudes are denoted by /ip with the sum over the index 
i running through all the possible bound states allowed by the Hamiltonian. These states 
are the color singlet qq elementary excitations identified in this work. The basic idea of NQA 
is to use the field equations of motion and derive integral equations for the Haag amplitudes 
and solve them to obatin a solution to the equation of motion. 

In order to solve for the Haag amplitudes it is necessary to calculate the mass and the 
coupling constant for each of the bound states. This has been accomplished successfully at 
finite T and ii for the two fiavored 't Hooft interaction model [0. In addition to bound state 
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masses and coupling constants it is also possible to determine the widths of these states. 
This quantity is the key to studying the confinement-deconfinement phase transition within 
our formalism. In the confined phase the bound states will have vanishing widths while in 
the deconfined phase we expect to see unbound resonant states with finite widths. Hence 
we propose to use the widths of the qq states as an order parameter to study the nature of 
the deconfinement phase transition within strong coupling QCD. 
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APPENDIX A: PROPERTIES OF FREE WILSON FERMIONS 

In this Appendix we present the properties of free Wilson fermions on the lattice in the 
Hamiltonian formulation [^. The free lattice Dirac Hamiltonian is given by 

X,l X 

~ J2 [^^(^ )7o*(^ + fii) + ^^(f + rii)7o*(x )] (Al) 

x,l 

where the third term is the Wilson term. For r = there is an eightfold fermion multiplicity 
which is removed when r 7^ 0. At each lattice site the free Dirac field in configuration space 
is given by 

p 

with u denoting the Dirac index. The excitation energy uj{p) will be determined shortly. 
The annihilation operators b and d annihilate the non-interacting vacuum state | 0). For 
our purpose it is not necessary to know the structure of the spinors ^ and rj. 

The only assumption that we shall make is that the creation and annihilation operators 
obey the free fermion anti-commutation relations 

[bHp), m] + = [d\p), d{q)] ^ = 5^,^ (A3) 
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Using this assumption we can recover the anti-commutation relations for the field operators 

[^,{t,x),^l{t,y)]^ = 6£,y-6,, (A4) 
provided that ^ and rj satisfy the relation 

Up)^lip) + Vp{-PH{-P) = S,. (A5) 
We normalize the spinors by demanding that the number density is given by 

AT = ^ : ¥{t,x)^{t,x) : = 2 ^ {b\p)b{p) - d\p)d{p)) (A6) 

X p 

where the symbol : : denotes normal ordering with respect to | ) and the factor of 2 
accounts for the spin degrees of freedom. Eq. (|A6|) fixes the normalizations of C, and rj to be 



UpKUp) = vAp)vlip)=2 (A7) 
Cip)Vui-p) = vlip)U-p) = (A8) 



which are consistent with Eq. 

In momentum space the charge conjugaton symmetric form of is 



P I 



-p] 



(A9) 



where the momentum dependent mass term is given by 

M(p) = M - r^cos(p- n;) (AlO) 
I 

The free Dirac field now becomes 

-^.it,p) = Kp)Up)e-'''^^^' + t/t(_p')^^(_p')e+-(?^)* (All) 
which is used to derive the equation of motion corresponding to Eq. (|A9|) 

ii>{tj) = : [vl>(t,p),i/0]_: (A12) 
= (^^sin(p-n07o7/ + M(p)7o^^(t,p) (A13) 

From Eq. ( |A13| ) one obtains the excitation energy 

/ \ 1/2 

= ( ^ sin2(p ■ ni) + M^{p) J (AM) 
25 



and the equations of motion for the ^ and 77 spinors 



UJ 



u{p) r]{-p) = -(^^sm{p-ni)-fo-fi + M{p)-fo^r]{-p) (A16) 

When r = these equations of motion are relativistic near the eight corners of the 
Brillouin zone denoted by ttq = (0,0,0), = (vr, 0,0), vfy = (0,7r, 0), vf^ = (0,0, vr), VTxy = 
(7r,7r, 0), TTxz = (7r,0,7r), TTy^ = (0,7r,7r) and TTxyz = (7r,7r,7r). The excitation energies near 
these values of momenta are equal which corresponds to the eightfold multiplicity mentioned 
above. This degeneracy is lifted when r 7^ due to the momentum dependent mass term 
Eq. (|A10|) . Using the equations of motion for ^ and 77 it is a simple excercise to show that 
the off-diagonal Hamiltonian vanishes and that the vacuum energy is given by 



{0\H^\0) = -2Vj2'^iP) 



(A17) 



Finally, we construct positive and negative energy projection operators A^{p) and A (p) 
as follows 



A+(p) = e(p)®e^(p) = ^ 



A (p) = r]{-p) ® rj\-p) = - 



1 + -j^ sin(p ■ n,)7o7i + 
1 ^ Vsin(p-ni)7 



M{p) 
uj{p) 



■7o 



M{p) 
'oil 7^70 



(Al8a) 
(Al8b) 



Note that the projection operators obey the condition 



[^\p) + ^-{p)\^, = 5,. 



(A19) 



as is required by Eq. (|A5|) 
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